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ABSTRACT 

We demonstrate that the mass of the most massive star in a cluster correlates non- 
trivially with the cluster mass. A simple algorithm according to which a cluster is 
filled up with stars that are chosen randomly from the standard IMF but sorted with 
increasing mass yields an excellent description of the observational data. Algorithms 
based on random sampling from the IMF without sorted adding are ruled out with a 
confidence larger than 0.9999. A physical explanation of this would be that a cluster 
forms by more-massive stars being consecutively added until the resulting feedback 
energy suffices to revert cloud contraction and stops further star formation. This has 
important implications for composite populations. For example, 10* clusters of mass 
10^ Mq will not produce the same IMF as one cluster with a mass of 10^ Mq. It also 
supports the notion that the integrated galaxial IMF (IGIMF) should be steeper than 
the stellar IMF and that it should vary with the star-formation rate of a galaxy. 

Key words: stars: formation - stars: luminosity function, mass function - galaxies: 
star clusters ~ galaxies: evolution - galaxies: stellar content - Galaxy: stellar content 



> 1 INTRODUCTION 



The insight that clustered star formation may be the domi- 
nant mode for star formation has grown over the last years. 
The form of the true distribution of stellar masses within 
these clusters, of the stellar initial mass function (IMF) , has 
been a subject of debate for a long time. The evolution of 
the stars, unresolved binaries, and the dynamical evolution 
of the clusters complicates the observational efforts to ex- 
tract the IMF. Unfortunately the most promising objects, 
very young stellar clusters (age < 3 Myr), are often still em- 
bedded in their natal cloud - again aggravating observations. 

Nevertheless the distribution of stars in young clusters 
seems to be fairly well described by a multi power-law func- 
tion with a slope or index (a) of 2.35 (the so-called ' Salpeter' 
value ) for stars with a mass larger than 0.5 Mq llKroupal 
I2OOII) . The initial mass function (IMF), 



useful to describe the stellar IMF with an invariant, 
multi-power law forrn (iKroupa et al.l 119931: iKroupa .2001 
iReid. Gizis fc Hawlevll20o3) ! 



ao = +0.30, 0.01 s; m/Mo < 0.08, 

ai = +1.30, 0.08 s; m/Mo < 0.50, 

a2 = +2.35, 0.50 m/M© < 1.00, 

as, = +2.35, 1.00 m/M©. 



(2) 
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where ^{m)dm is the number of stars in the mass in- 
terval m, m + dm. Several observations find the Salpeter 
value (as = 2.35) for a large variety of condl 



tions ijMa ssev 
Parker et al.i i2 
Bell et al.l I2OO 



Huntei' 1998; 



riety 

Siri anni et fai] I2OOL 
Iwvsee^l] 



Parker et_al.ii200k iMassevi i200a. 

iPiskunov et al.l 



2002 



2002 



It is therefore 



e-mail: cweidner/pavel@astro.uni-bonn.de 



We refer to this form as the standard or canonical stel- 
lar IMF because this form fits the luminosity function of 
Galactic-field and cluster stars below 1 M© and also repre- 
sents young popula tions above 1 Mq iKroupa et al.lll993l : 
lKroupall200ll . l2002l) . As pointed out bv lScald ^19'9al2005^ 
though, significant uncertainties remain in the determina- 
tion of the IMF to the point that the case can also be 
made that a single form of the IMF may not exist. In view 
of this, the ansatz made here and elsewhere is to propose 
the hypothesis of an invariant standard or canonical IMF 
(eq. 121 and to test if the variation of the observed IMF can 
be understood to be the result of astrophysical effects (ob- 
scuration, stellar evolution, stellar multiplicity), dynamical 
effects (mass segregation, stellar evaporation and ejections), 
stochastic effects (finite A'^-sampling from the IMF) and the 
construction of composite populations (addition of many dif- 
ferent clusters). 

Similarly, the embedded cluster mass function (ECMF) 
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has been found to be weU-described by at least one power- 
law, 

^ocl(Mecl) OC M-f , (3) 

where dA'"cci = £,cc\[Mcc\) dMcd is the number of 
embedded clusters in the mass interval Meci, Med + 
dMeci and Med is the cluster mass in stars. The 
observational evidenc e points to a pos sibly universal 
form of the ECMF: iLada fc Ladal tOO^ find a slope 
/3 = 2 in the solar neighbourhood for clusters with 
masses between 50 and 1000 Mq, while iHunter et alJ 
l|2003h find 2</3<2.4 for 10^ < M e d /Mg < lO'^ in the 
SMC and LMC, and IZhang fc Falll l)l999l^ find 1.95 ± 
0.03 for 10*^ <Med/Mf7i < l O'^ in the Antennae galaxies. 
I Weidner. Kroupa fc LarsenI ll2004ll discovered that /3 = 2.35 
best reproduces the observed correlation between the bright- 
est young cluster and the galaxy-wide star-formation-rate 
for a large sample of late-type galaxies. 

As already mentioned by Vanbeveren' and dis- 

cussed in more detail by Kroupa & Weidner ( 20q3), the 
composite or integrated galaxial stellar IMF (IGIMF) is ob- 
tained by summing up the stellar IMFs contributed by all 
the star clusters that formed over the age of a galaxy, 

/•A4ol,max(SFfl(t)) 
ClGIMF(m;f) = / ^ "imax(Mcd)) 

•^-M"<!cl,min 

■eod(Med)rfMed, (4) 

where Ccd(Med) is the ECMF and ^(m < (Med)) is 

the stellar IMF in a particular cluster within which the 
maximal mass of a star is mmax. Med,min {~ 5Mq, Taurus- 
Auriga-type "clusters") is the minimal cluster mass, while 
the maximal cluster mass, Med. max, depends on the galaxy- 
wide star-formation rate, SFR iWeidner. Kroupa fc LarsenI 
l2004l) . 

A critical function entering this description is thus 
fWinax(Med). Assumiug the stellar IMF to be a continuous 
distribution function, this mass of the most massive star in 
an embedded cluster with the total mass Med in stars is 
given by 

1 = / ^{m)dm, (5) 

with 

Mod = / mS_{m)dm, (6) 

since there exists exactly one most massive star in each 
cluster, and neglecting statistical variations. Here miow 
— 0.01 Mq is the minimal fragmentation mass and 
ntmax* ~ 150 M|T) is the measured maximal stellar mass 
limit (|Weidner fc Kroupal 12004: Fieeil 1200,4 lOev fc Clark3 
I2OO5I) . On combining eqs. |^and|H|the function 

mmax = 'mi^ax(Med) (7) 

is quantified bv lWeidner fc Kroupal l)2004h and in §|21 This is 
the analytical ("ana") maximum-stellar-mass-cluster-mass 
relation which incorporates the fundamental stellar up- 
per mass limit (noted by the leading superscript "1") of 
m,max« = 150 A/0. Later-on other maximum-stellar- mass- 
cluster-mass relations are indicated by different superscripts: 



"ran", "con" and "sort" for the different Monte-Carlo sam- 
pling methods (see 12.211 and also "u" for the case without a 
fundamental stellar upper mass limit. 

Weidner fc Kroupa (20o3) infer that a fundamental 
upper stellar mass limit, rrimax* ~ 150 Mq, appears to 
exist above which stars do not occur, unless 03 ft 2.8, 
in which case no conclusions can be drawn based on 
the expected number of massive stars. As reviewed by 
iKroupa fc WeidneJ ((200^, the existence of such a stellar 
upper mas s limit has been furt her substantiated by iFigeil 
J2005I) and lOev fc Clark j (|200^ for a range of star clusters 
and OB associations. 

We thus have, for each Med, the maximal stellar mass, 
mmax(Med) ^ JTimax* , aud with this in formation eg . IH can 
be ev aluated to compute the IGIMF. iKroupa fc Weidneil 
J2003h find the IGIMF, when evaluated to the highest 
cluster masse s, to be significantl y steep er than the stel- 
lar IMF, and I Weidner fc Kroupal (l2005ti extend the anal- 
ysis to a time varying ECMF by noting that Med, max in- 
creases with the star-formation rate of a galaxy. They show 
the IGIMF to be not only steeper than the stellar IMF, 
but also to depend on galaxy type. The implications of 
these findings are rather significant for the supernova rate 
jWeidne r <^ Kr9x^p^|200, 4) and for the chemical evolution of 
galaxies iR^ppen. W eidner fc Kroupal l2006l . 

But these r esults remain not without a challenge. 
lElmegreenI l)2004h argues that there is no evidence of a re- 
lation mmax = mmax(Mcd) ^ mmax*. This relation im- 
plies that many small, low-mass, star-forming events will not 
have the same combined IMF a s one major SF even t of th e 
same mass. Thus, according to iKroupa fc Weidneil i2005l) . 
10^ clusters each with a mass of 20 Mq would provide a 
combined IMF that differs from that of one cluster with a 
mass of 2 X 10® Mq by being underrepresented in stars with 
a mass above about mmax = 1 Mq . The contrary, often 
voiced view is that stellar masses sample the IMF purely 
statistically such that 10^ clusters containing 50 stars (on 
average 20 Mq) will give the s ame combined IMF as one 
cluster containing 5 X 10*^ stars llElmegreenlll99gl. 120041) . 

With this contribution we demonstrate conclusively 
that the purely statistical notion is false, and that the 
stellar IMF is sampled to a maximum stellar mass that 
correlates with the clust er mass. Therewit h we affirm 
the results obtained by iKroupa fc Weidneil i2003ll and 
[Woidncr & Kroupa ( 20o|), and we also attain useful insights 
into the process of star-cluster formation. 

In Section (21 our Monte-Carlo procedure is described 
and the maximal-star-mass-cluster-mass relation is derived, 
while in § 13 the Monte-Carlo experiment is applied to the 
IGIMF and the results are presented. A discussion with con- 
clusions is available in §|31 

2 THE MAXIMAL STAR MASS IN A 
CLUSTER 

2.1 Previous' studies 

Over the past 20 years several studies investigated a possible 
connection of the maximum stellar mass in a cluster and 
the mass of the cluster because such a relation, if it were to 
exist, would allow important insights into the star-formation 
process. 
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iLarsorJ (^S^) compared the properties of molecular 
clouds with the spatial distribution of the associated stellar 
population. He found that more massive and dense clouds 
favour the formation of massive stars and fitted the following 
formula to the observations, 



= 0.33M°i«d- 



(8) 

He re-evaluated jLarsonl 12063) this equation with more re- 
cent data and applied it directly to the cluster mass instead 
of the cloud mass, 

mn,ax = 1.2M°■f^ (9) 

Thi s correlation is plo tted as a dotted line in Figs. and Q 
lElmegreenI il983l) proposed a model for the origin of 
bound galactic cluster s wher e the luminosity of the stars 
from a iMiller fc Scald Jl979D IMF overcomes the binding 
energy of a molecular cloud core. The star-formation effi- 
ciency then discriminates between bound clusters and OB 
associations. He derived an analytical estimate for a relation 
between the maximal star mass and the cluster mass from 
statistical considerations regarding the appearance of stars 
with various masses from the Miller-Scalo-IMF, 



Mod = 



e-*2F[yli(logm„ 



A. 



1 -F[Ai(logmn 



C2)y 



with 



(10) 



(11) 



and Ai = {2CiY'\ ^2 = In 10(C2 -hlnlO/4Ci) = -1.064, 
A'i = C2+ lnl0/2Ci = 0.065, Ci = 1.09 and C2 = -0.99. 
In Figs.0and|7|this relation is shown as a short-dashed line. 

On the othe r hand, using a single po wer-law IMF with 
a lSalpeteil ^1955^ slope. lElmeereenl ll200Cf) solved eqs.|Sland 
|S| assuming mmax* = 00, with the result. 



Mod = 3 X 10"^ 



lOOM, 



(12) 



whi ch is shown as a long-dashed line in Figs. U and 13 

iBonnell. Bate fc Vhi3 ll2003f) and lBonnell Vine fc Bat3 
l|2004il numericallv studied star-formation in clusters us- 
ing their smooth-particle-hydrodynamics (SPH) code. Here 
a turbulent molecular cloud fragments hierarchically into 
smaller subunits. When the density of a clump gets higher 
than a critical value, it is replaced by a so-called 'sink' par- 
ticle which only lets matter in but not out. These sink par- 
ticles form the final stellar cluster by interactions and merg- 
ings and are called stars at the end of the simulation. This 
hierarchical cluster formation scenario leads to the relation, 

mmax OC , (13) 

and is shown as a short- dashed- dotted line in Figs. and 
[7| There eq. 1131 is normalis ed to nzmax = 27 for a 
cluster of Med = 580 M r iBonnell. Bate fc Vind [20031: 
iBonnell. Vine fc Batj|2004f) . It should be noted here that 
these simulations do not include magnetic fields, stellar feed- 
back and stellar mergers, all of which are believed to be of 
gre at importance for star-form at ion . 

IWeidner fc Krou pal (l2004^ started with similar assump- 
tions"as ^hnB^^enril200ol) but included a physical upper 
limit for the stellar mass, rrimax* = 150 M©, while solving 




Figure 1. The maximal star mass, mmax, in depend ence of 
the cl uster mass, Mgd, both in logarithmic units. The iLarsonl 
('2003^-rolation (eq. ^ is shown as a dotted line, while the 
Elmcercen 1 1983)-rclation (oq. llOi is the short-dashed line. 
The Elmegreeni 1 2000)-rclation (eg. 1121 is the long-dashed line. 
The resul t (eq. [131 from numerical SPH star-formation sim- 
ulations iBonnell. Vine fc Batd l2004h is plotted as a short- 
dashed- dotted line an d as a triangle for a specific model from 
IBonnell. Bate fc Vind l2003l) . while the long -dashe d- short- dashed 
line marks the expectation values (eq. 1141 from lOev fc Clarkel 
12005) an d the thick solid line sho ws the semi-analytical model 
(eq.0 of iWeidner fc Kroupal <2004) for the standard IMF (eq.|2j. 
The thin solid line shows the identity relation, were a "cluster" 
consists only of one star. 



eqs. 13 and |S| Due to mmax* and using the standard multi- 
power law IMF (eq.|21l, the result cannot be written out an- 
alytically but the equations have to be solved numerically. 
The result is plotted as a thick solid l ine in Figs. U and [ 
In a broader statistical analysis, lOev fc Clarkel (|2(1 
calculated the probabilities that the observed upper mass 
limits in several clusters and OB associations in the MW, 
LMC and SMC come from a sample with a fundamental 
upper mass limit, mmax*, or not. They conclude that an 
upper mass limit between 120 and 200 Mq is the most likely 
explanation for the observed maximum masses. In oder to do 
so, they calculated the expectation value for the maximum 
mass if a number A*' of stars is randomly sampled from an 
IMF. This yields the equation. 



(mmax) = mn 



^(m)dm 



dAfod. (14) 



Integrating this numerically yields the long-dash-short- 
dashed line in Figs. and |7| assuming mmax* ~ 150 M©. 



2.2 Monte-Carlo experiments 

All the above mentioned investigations suggest that the clus- 
ter mass indeed appears to have a limiting influence on the 
stellar masses within it. However, it would be undisputed 
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that a stochastic sampling eflect from the IMF must be 
present when stars form. To investigate the possible exis- 
tence of a maximal stellar mass in clusters statistically and 
to confirm or rule out if such a relation is purely the result 
of the random selection from an IMF, three Monte-Carlo 
experiments are conducted: 

- pure random sampling (random sampling), 

- mass-constrained random sampling (constrained sam- 
pling), 

- mass-constrained random sampling with sorted adding 
(sorted sampling). 

For each, two possibilities are probed: stars are sampled from 
the IMF without a maximal mass (mmax* = cxd^), or their 
masses are limited by mmax* = 150 Mq . 

2.2.1 Random sampling 

For the random sampling experiment, 2.5 x 10^ clusters are 
randomly taken from a cluster distribution with a power- 
law index /3jv = 2.35. The clusters contain between 12 stars 
and 2.7 x lO'^ stars. The relevant distribution function is the 
embedded-cluster star- number function (ECSNF), 

diVe.l (X iVf,~ „ (15) 

which is the number of clusters containing N £ [A'^' , ' -|- 
dN'] stars. Each cluster is then filled with A*' stars randomly 
from the standard IMF (eq. |5J without a mass limit, or 
by imposing the physical stellar mass limit, m ^ 150 M©. 
The stellar masses are then summed to get the cluster mass, 
Meci. Note that such a cluster distribution gives an embed- 
ded cluster mass function (ECMF) that is virtually identi- 
cal to eq.|^(Fig.|21l. The resulting distribution of maximum 
stellar masses is plotted in the mmax, M^d plane (Fig.|3J as 
contour lines, to show the overall distribution (for a more 
detailed discussion see S 12. 2. 41 . 

As in this method the higher cluster masses are only 
very scarcely sampled, a second method is used to evaluate 
the mean of maximal masses in more detail. In order to 
do so the cluster star numbers A'' = 12 to 2.7 x 10^ are 
divided into 10 logarithmically equally-spaced values. Each 
of these numbers is then filled 10000 times with stars form 
the IMF, keeping only the mass of the most massive star for 
each cluster. The mean maximum mass for every bin then 
defines, 

^^max (A/ecl), (16) 

with a limit of mmax* = 150 Mq and 

"'Wmax(A/ccl), (17) 

in the unlimited case. 



2.2.2 Constrained sampling 

In this case 2.5 x lO'^ clusters are randomly taken from 
the ECMF (eq. 01 between 5 Mq (the minimal, Taurus- 
Auriga-type, star-forming "cluster" counting ~ 15 stars) 

^ For practical reasons of numerical computation rrtmax* is 
adopted to be rrimax* = 10^ Mq in the unlimited case. 




log,„(M,yMj 



Figure 2. The embedded cluster mass function (ECMF) derived 
from randomly picking lO'^ clusters using ec\.^^(dashed line) and 
as described in ^ 12.2.11 (solid line). The slopes are virtually the 
same (/3jv = /3)- Only for very small cluster masses does the 
solid line deviate due to the underlying standard IMF because 
small- A'^ clusters can nevertheless have masses > 10 Mq if their 
constituent stars happen to be massive. This accounts for the turn 
down below 10 Mq and the surplus in the range 10-150 Mq. 

and 10^ Mq (an approximate maximum mass for a sin- 
gle stellar populati on that c onsists of one meta llicity and 
age, Wei dner. Kroupa fc Larse n 200^; ICieles et al. 2Q05i) 
and again with /3 = 2.35. Note that /?jv ~ P because the 
ECSNF and the ECMF only differ by a nearly-constant av- 
erage stellar mass (Fig. |5J . Then stars are taken randomly 
from the standard IMF and added until they reach or just 
surpass the respective cluster mass, Med. Afterwards the 
clusters are searched for their maximum stellar mass (plot- 
ted as contours in Fig. |31 see S l2.2.4l for discussion). 

Again the sampling of high cluster masses is very poor. 
Therefore, the cluster masses from 5 Mq to lO*' Mq axe di- 
vided into 10 logarithmically equally-spaced values. Then 
each of these 10 cluster masses is filled 10000 times with 
stars form the IMF until their combined mass reach or just 
surpass the cluster mass and only the maximum star mass is 
recorded. The average mmax values for each of the 10 cluster 
masses define the relation, 

'mrx(Med), (18) 
in the limited case and 

"mrx(Me.i), (19) 
in the unlimited case. 

2.2.3 Sorted sampling 

For the third method again 2.5 x lO'^ clusters are randomly 
taken from the ECMF (eq. between 5 Mq and 10** Mq 
and with /3 = 2.35. However, this time the number A'' of 
stars which are to populate the cluster is estimated from 
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A'^ = Meci / m-av, where rriav = 0.36 is the average 
stellar mass for the standard IMF (eq.|5J between 0.01 Mq 
and 150 Mq. These stars are added to give Mod, ran, 

Mod, ran — mi, 
N 

such that rai mi+i. If Mran < Mod in this first step, an 
additional number of stars, AN, is picked randomly from 
the IMF, where AN = (Mod — Med, ran) / n^av Again these 
stars are added to obtain an improved estimate of the desired 
cluster mass, 



*Mo, 



N+AN 



again with rrii ^ m^+i. When *A/od,ran surpasses A/od, it is 
checked whether the sum is closer to Mod when the last star 
(the most massive one) is subtracted or not. If * Med, ran is 
closer to Med with the last star this one is the most massive 
one for the cluster, otherwise it is the second last star ^ . This 
procedure is followed and repeated until all clusters from the 
ECMF are 'filled'. The contour plots of the most massive 
star for each cluster are shown in Fig. |3 and discussed in 
§11111 

Again, for a more detailed analysis, 10 cluster masses 
are generated like in § l2.2.2l but filled with stars in the sorted 
way described above. The mean over every of the 10 cluster 
masses then yields the relation 



I sort 



(Mo. 



(20) 



2.2-4 Comparison of the samplings 

In order to exemplify the difi'erences between the three 
sampling methods the following gedanken experiment may 
be considered: 

A sample of 10 stars consists of 9 stars with 5 Mq and one 
with 11 Mq. For the random sampling, this would be a 

cluster with Med, ran ~ 56 Mq, with mmax.ran ~ H Mq. 

If, for the sorted sampling, the aimed-at cluster mass is 50 
Mq, the actual cluster mass would be 45 Mq {= 9x5), 
because 45 is closer to 50 than 56, and mmax,sort would be 
5 Mq. In the case of constrained sampling the order would 
be important. If the aimed-at cluster mass is 50 Mq and 
the 11 Mq star is among the first 9 stars, Afod.con would 
be 51 (=8x5 + 11) and mmax.con = 11 Mq. But if the 11 
Mq star is the tenth star, Med, con would be 45 Mq and 
^Tlniax, con — 5 Mq as in the sorted-sampling case because, 
45, rather than 56, is closer to 50 Mq. 

Figs. 13 11| and 13 plot the contour lines of the most mas- 
sive stars of all the simulated clusters for random sampling 
(Fig. 13 , constrained sampling (Fig. 0J and sorted sampling 
(Fig. 13, all with the physical limit mmax* = 150 M0. 

As can be seen from Fig. |3 the random sample occu- 
pies the whole parameter space between the two extremes, 
which are either that nearly the whole cluster consists only 



^ For example, if Afed = 50 Af©, mjv 
5Mq, then for *A/oci,ran = 52 Mq, m 
52 is closer to 50 than 42. If Afgd would be 45 Mq, mmax would 
be 5 Mq. 



lOA/0 and m^f—i 
max would be 10 Mq, 



1.5 



o 




0.5 



1 2 3 4 5 6 

logic (Meci/M,,) 

Figure 3. Maximal mass of stars versus cluster mass. The 
contour lines show how often a certain combination of clus- 
ter mass and maximal star mass occurs. The further-out they 
lie (towards the upper right) the smaller is the number of 
clusters with this combination. They are spaced logarithmi- 
cally such that the outer-most one represents a single clus- 
ter with a certain mass and maximal star (logj^Q A^ocl = 0), 
eg. Afed = 10^ Mq, mmax ~ 100 A/q. The inner-most one 
(near A/ed = 7 Mq, mmax ~ fewM0) stands for log^^Q A^ecl — 
4.5 clusters with this maximal star mass. The lines in-between 
are in steps of 0.5 dex. Here the contour lines are the result of 
the Monte-Carlo simulations with random sampling. The semi- 
analytic result (eq. [7J, 'm^^^{Mc^i), is the thick solid line. Mean 
values are shown as the thick dash- dotted line for random sam- 
pling ('mmax(A^ocl), ^1- lltil . as the thin dashed line for con- 
strained sampling ('mmax(A/ed), eg. 1181 and as the thin dotted 
line for sorted sampling ('m^ax(A^ecl), eg. 1201 . The identity re- 
lation mmax = A/ed is plotted as a thin solid line. 



of low-mass stars, or a single star accounts for the entire 
cluster mass (mmax = M^d). Such one-star-clusters would 
correspond to freak star-formation events such as is envis- 
aged for_ava£iable_gas;^ec^^ the- 
ory ihi. Mac Low fc Klessenll2005^ . For sorted sampling the 
clusters are shifted more towards smaller stellar masses and 
never touch this line. While the constrained sampling lies in 
between these two extremes and still populates the param- 
eter space up to the identity relation. 

Note that 'mmix(Med), shown in Figs.|3to|3(f/iicfc dot- 
ted line in Fig.|3 thin in Figs. |3 and |1J, is nearly identical to 
the semi-analytical estimate 'rn^^^{Mcci) (eq. thick solid 
Ime in Figs.|3to|3. The slight deviations are probably due 
to the stochastic element in the process to decide which is 
the most massive star in the sorted Monte-Carlo experiment 
(see S iTOl . 

The non-smoothness of the contour lines in the upper 
right edge of the figures is an effect of low-number sampling. 
Here only the outer-most contour line is present - indicating 
that there are only single events in this region of parameter 
space obtained with 2.5 x 10^ clusters, which have together 
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1.5 




0.5 



1 2 3 4 5 6 

logio (Meci/M^) 

Figure 4. Like Fig.|3]but this time the contour lines are the result 
of the Monte-Carlo simulations with mass-constrained sampling, 
and 'mJj!['"^(Moci) is shown as a thick line, while 'r?i™jj(Meci) and 
'"*max(-'^ccl) are thin, 'rrajjjljj is the thick solid line. 



1.5 




0.5 



1 2 3 4 5 6 

logio (Mec/Mg) 

Figure 5. Like Fig. 1!^ and Bl but this time the contour lines are 
the result of the Monte-Carlo simulations with sorted sampling 



about 600 • lO'' stars. This number of clusters comprises the 
computational limit of the available hardware. 

The agreement of the m ean value 'm^^yfA^ ec i) an d 
the semi-analytic result of IWeidner fc Kroupal i2004l . 
'"^max(Afcci), is iu principal not surprising. The method of 
sorted adding of stellar masses in order to get the cluster 
mass corresponds to a Monte-Carlo integration of eqs.j^and 



1^ Therefore the result should agree with the numerically in- 
tegrated (semi-analytical) one. 

Another difference between the samplings lies in the 
average mean stellar masses in the clusters. To determine 
these, cluster masses from 5 Mq to 10*' Mq are divided in 
10 logarithmically equally-spaced values and each is filled 
10000 times with stars from the IMF in the three different 
ways described in S 12.2.11 r2.2.2l and l2.2.3l The mean stellar 
mass for each cluster is calculated by 



JV 



(21) 



where A'^ is the number of stars in each cluster. Then for 
each cluster mass, Mcdj (j ~ 1...10), the 'mean of means' 
is computed by 



m-j(A'/oci,j) = in— 



1 

Nccl 



JVnH 



(22) 



with iVcci = 10000. 

The different average mean stellar masses are shown in 
Fig. |S| For the random sampling the mean is always about 
0.36 Mq, as expected for the canonical IMF between 0.01 
and 150 Mq . The other sampling methods have lower means 
for light clusters which rise up to 0.36 Mq for more massive 
ones. This is a result of the limit which low-mass clusters im- 
pose on their stellar content. For random sampling only the 
number of stars determines the cluster mass, in contrast to 
reality where the natal cloud-mass and the star-formation- 
efficiency rule the final cluster mass. 

Also shown in Fig. |S| is the average mean stellar mass 
for the analytical model (eq. |7J , which is given by 



™max(A/ecl) 



m(^dm 



where ^(m) is eq. |5| 

2.2.5 Comparison with observational data 

Table □ shows a non-exhaustive compilation of cluster 
masses and upper stellar masses for a number of MW and 
LMC clusters (see appendix ^ for more details). These 
data are shown as dots with error-bars in Figs. |7| |S| 
and El The result of a clu ster formation simulation by 
iBonnell. Bate fc Vind ll2003h is shown as a large triangle. 

Fig. |7| compares the mean mmax values of the Monte- 
Carlo experiments together with the observations from 
Tab. [Tl the semi-analytical result of IWeidner fc Kroupal 
I2OO4I) and the different results of the previous' studies 
(S I2.III . The "unlimited" mean values are marked with a "u", 
while the limited ones (mmax* ~ 150 Mq) with an "1". The 
mean values of the "u" cases are all clearly distinct from the 
observations and not regarded further, while for the "1" cases 
the mean values are in reasonable agreement with the obser- 



(23) 



and ''"mS?^x(M,ci) lie 



vations. But especially ''''m™x(Meci 
above the observations, while 'mmix(-Mcci) and 'mmax(M3ci) 
fit the observations rather well. 

The ob servational analysis (eq. |H| and O b^jLaraoi] 
l|^l|2003) and the analytical result (eq.jTSl of lElmeeree" 
C19S3) have a shallower slope and underestimate the ob- 
served stellar masses for Med > 100 Mq . The analyt- 
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0.35 - 



.0.3 



0.25 



0.2 




Table 1. Empirical cluster masses, maximal star masses within 
these clusters and cluster ages from the literature. 



3 4 
log,„(M^„,/Mj 

Figure 6. Average mean stellar mass for the different Monte- 
Carlo samples against cluster mass (eq. 1221 . For random sam- 
pling (dotted line) the value is constant within the numerical 
noise around the expected value of 0.36 Mq. In the case of con- 
strained sampling {dashed- dotted line) and sorted sampling {long- 
dashed line) it starts rather low for low mass clusters and rises 
up to expected value of 0.36 Mq. The solid line shows the rela- 
tion described by eg. 1231 Note that the long- dashed line lies above 
the solid curve because the maximal stellar mass in the sorted- 
sampling algorithm is systematically higher than the analytical 
result (Fig. |7| to |5] below) for cluster masses below 10* Mq. 



ical result (eq. I12|l of 'ElmeEreen' ('20 Q(]|) and the star- 
form ation si mulation (eq. | 13[l by IBonnell. Bate fc Vinel 
EmZ ): Bon nell. Vine fc Bate! (l2004) fit much better but ig- 
nore the obse rved upper mass limit f or stars n ear 150 Mq 
JWeidner fc Kroupa 2004: .Figer. 200 j : lOev fc Cl arke 20^ . 
TaJien this int o account lead s to th e semi-analytical (eq. |7|l 
relation. The iQev fc Claris ||200^ expectation values for 
TUmax (eq. I14II follow quite well the constrained sampling re- 
sult, but in doing so they over-predict mmax in comparison 
with the observations. 

In Fig.|51only "rn;^ax(Ai"cci) and 'mmax(Afcci) are shown 
but for different slopes of the IMF above 1 Mq . For 
"^max(A/oci) 03 = 2.35, 2.70 and 3.00 is plotted, while 
for ' rri^ax( A/eci) the relations are plotted for 0:3 = 2.35 
and 2.70. In the case of "m5j!J'"x(Afeci) we can see that the 
steeper IMF slopes fit the obser vations much better. There- 
fore, if ISagar fc Richtlerl il99l[) are right about their no- 
tion that unresolved binaries would mask a steep underlying 
IMF, constrained or even random samplin g cannot be ruled 
out. B ut it should be noted here, that the lSagar fc RichtleJ 
il99ll) examination was only carried out for stars with 
masses between 2 and 14 Mq. Therefore, the O-star regime 
has not yet been explored in this respect. 



2.3 Ageing of the stars 

To see if stellar evolution together with constrained sam- 
pling can mimic the effect of the sorted sampling, stars with 



Designation 




^maxobs 


age 


Source 






\ Ti/T 1 

We\ 


[Myr] 




Tau-Aur 


25 ± 15 


2.2 ± 0.2 


1-2 


(1) 


Ser SVS2 


30 ± 15 


2.2 ± 0.2 


2 


(2) 


NGC1333 


80 ± 30 


5 ± 1.0 


1-3 


(3) 


p wpii 


1 on -t- 


s -t- 1 n 


U. 1- J- 




IC348 


109 ± 20 


6 ± 1.0 


1.3 


(5) 


NGC2024 


225 ± 30 


20 ± 4 


0.5 


(6) 


a Ori 


225 ± 30 


20 ± 4 


2.5 


(6) 


Mon R2 


259 ± 60 


10 ± 1 


0-3 


(7) 


NGC2264 


355 ± 50 


25 ± 5 


3.1 


(8) 


NGC6530 


815 ± 115 


20 ± 4 


2.3 


(9) 


Ber 86 


1500 ± 500 


40 ± 8 


2-3 


(10) 


M42 


2200 ± 300 


45 ± 5 


<1 


(11) 


NGC2244 


6240 ± 124 


70 ± 14 


1.9 


(12) 


NGC6611 


2- 10* ± 10000 


85 ± 15 


1.3 ± 0.3 


(13) 


Tr 14/16 


4.3-10* -2.3-10*/+2-10* 


120 ± 15 


<3 


(14) 


Arches 


540* -3.5-10V+2-10'' 


135 ± 15 


2.5 


(15) 


R136 


140^ -5-10*/+1.5-10S 


145 ± 10 


1-2 


(16) 


Simulation 


580 


27 


0.456 


(17) 



^iBricefio et alj <2002l) ^iKaas 1 



Asgii 



<2003ft: 



.^...^^ et"ai] i2004^. 

ICetman et alj ^2002^ . * IWilking Lada fc Younj <1989h: iLarsonl 
|2003^. '^iPreibisch fc Zinneckej <200ll):lLada fc Ladal i2003l) . ^ 



^herrv. Walter fc Wolld |200. 



JSung^ 



J? r2 



'^ ICarpenter et al.l il99l 

9 'Prisinzano et al.l |l20Q5 



B esscll & Chun' f2004') 

iDa miani et al. (2004), ^" Masscv. Johnson & Dc Gioia-Ea a^ood 
Eii^nyallcnari ct al. (1999), " Hillenbran d fc Hartmam] 

11998); Hillenbrand ct al. (1998), 

lMassey^~ ohnson & DcGioia-Eastwood 11995): Pa rk fc Sund 
^0021 . 13 Bonatto. Santos Jr fc Bica l.20Q5i). 1* 

iMassev fc JohnsonI <1993t): 

iMassev. Johns on^^DeGioig^asWoodl ^1995^. ^^iFige r et alj 
J2002^ I'^ lMassev fc Hunteil <1998l):ISelman et alj 119991) . 
iBonnell. Bate fc Vind ^2003^ 



m 50 M0 are ev olved with the single stellar evolution 
(SSE) package from iHurlev. Pols fc Toiiil (|S)0(1), while for 
more-massive stars fitting fo rmulae derived from s tellar evo- 
lution models computed by iS challer et alJ lll992l) are used 
(see Appendix IHl for for the detailed fitting functions). For 
this purpose 1 x 10^ clusters are chosen from an ECMF 
with (5 = 2.35 and then aged for 1, 2, 2.5, 3 and 3.5 Myr. 
The evolved stellar masses are added after excluding neu- 
tron stars and black holes to give Med and then searched 
for the most massive star in each cluster. 

In Fig. 1^1 the effect of this ageing is shown. Within the 
first 2.5 Myr the mean values of the sorted sampling and the 
constrained sampling algorithms are clearly distinct. The 
lines change due to mass loss of the heavy stars, which 
amounts to about 40 to 50 % of the initial stellar mass, 
according to the models used (see Appendix [CI for a com- 
parison of different models). After 2.5 Myr the massive stars 
start to explode as supernovae and the constrained sampled 
clusters (''rn'^^^{Mcc\)) move closer to the sorted-sampled 
ones ('rfCax(Meci)) for Meci > 10^ Mq. After 1 Myr and 
for clusters with Meci ^ 10^ Mq, ageing shifts 'mJj5'"x(Mcci) 
closer to the observations, making a distinction between con- 
strained and sorted sampling not possible. Nevertheless, the 
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1 2 3 4 5 6 



log,„(M,yMj 

Figure 7. The thick solid line shows the dependence of the 
maximal star mass on the cluster mass for as = 2.35 from 
the semi-analytical model ('mj^jj^(Moci), eq. |7J. The thick short- 
dashed line shows the mean maximum stellar mass for sorted sam- 
pling ('m^^{Mf,ci) , see S I2.2.iji . The long-dashed lines are mass- 
constrained random-sampling {^'^rn'^^^(M^^i), see i) l2.2.2l results 
with an upper mass limit of 10® Mq {straight line) and 150 Mq 
{curved line). Pure random sampling models (''"m™x(Meci), 
see § 12.2. li arc printed as dot- dashed lines. The curved one is 
sampled to mmax = 150 Mq while the straight one assumes 
'"max* = 10® Mq . The thin solid line shows the identity re- 
lation, were a "cluster" consists only of one star. The dots with 
error bars are observed clusters (see Tab.0, while the triangle is 
a res ult from a star-formation sim ulation (eq. Iliji with an SPH 
code iBonnell. Bate fc Vinell2003l). A i previous study of massive 
stars in clusters (ea.l51andl51 bv lLarso ij(1982, 2003j) is shown as a 
thin d otted line, while an analytic estimate ('ca. llOI _bv|ElmegreCTi| 
^1983^ as a thin short-dashed line and eQ. ll2l bv El rnegreein200cF 
as a thin long- dashed line. The thin long- d ashed- short- dashed lin e 
marks the expectation values (eg. 1141 from lOev &: Clarkj l2005l) . 



observational data for Med < 10^ Mq agree much better 
with sorted than with constrained sampling for all ages. 

It must be noted here that the observed stellar masses 
(Wmaxobs) in Tab. Q are a mixture of present-day (PD) 
masses and zero-age main-sequence (ZAMS) masses. For 
stars below roughly 50 Mq this is not critical, as for them 
mass-loss is not dominant during the first 3 Myr. But the 
three most-massive clusters have ZAMS maximal stellar 
masses and can therefore not be compared with the aged 
tracks in Fig.|Ul Our previous f? l2.2ll comparison of the data 
with zero-age main-sequence isochrones is thus justified. 



2.4 Statistical analysis 

In order to evaluate the statistical significance of the dif- 
ferences between the Monte-Carlo simulations and the ob- 
servations, two statistical tests are applied. The statistical 
tests are only applied to the zero-age Monte-Carlo samples, 
as the PD and ZAMS masses for low-mass stars do not dif- 



= 2.35 = 3.00 




1 2 3 4 5 6 

Figure 8. The thick solid line shows 'mj^jf^{Moci) for 03 = 2.35 
(eq. 0. The thick dashed lines are 'rn^Jx(^ecl) for 03 = 2.35 
(upper line) and 2.70 (lower line). The dot- dashed lines are con- 
strained sampling Monte-Carlo results, "mJjJ°^(Meci), with three 
different slopes {013 = 2.35, 2.70 and 3.00 above I Mq) of the 
input stellar IMF. The thin solid line shows the identity relation. 
The dots with error bars are observed clusters (Tab.lTl. while the 
triangle is a result from the star-form ation simulation with an 
SPH code iBonnell. Bate &: Vinell2003l) . 



2 



s 



' ' y ' ' '0 Myr 


' / 


1 Myr 

" " " 2.5 Myr 

3 Myr_ _ | 


^ / .4^^' 


3.5 Myr 






\ / \ 
- 1 '''fx 

.' '/ 
1 */ 




- ff^ 

1/ 

1 






1 2 3 4 5 6 

log,„(M,yMj 

Figure 9. As Fig. 171 but the mean curves include ageing by 1, 2, 
2.5, 3 and 3.5 Myr. The stars in the Monte-Carlo-simulations 
ar e subject to stellar evolutio n according to the SSE package 
bv lHurlev. Pols fc Toull i2000h and our own extensions for stars 
^ 50 Mq which includes not only finite life-times but also stellar 
mass-loss. The thick dashed lines are for clusters which are con- 
structed using sorted sampling, while the dot-dashed lines are for 
constrained sampling. Note that 'm^°3j(Afeci) = '"^axC^ccl) for 
ages ^ 3 Myr and Afod ^ 10^ Mq . 



Maximum-stellar-mass-cluster-mass-relation 9 



fer substantially for such young clusters and all the massive 
stars (m ^ 50A/q) in the observational sample are always 
ZAMS masses. While for ages ^ 2.5 Myr different stellar 
models agree rather reasonably on the properties and evo- 
lution of massive stars (see Appendix |n| for a substantial 
comparison of some models) the theoretical foundation of 
the formation of such massive stars is s till very weak (for a 
detailed discussion of massive stars, see lKroupa fc Weidneil 

The probabilities that the observed masses of the most- 
massive stars, mmaxobs, in the clusters are drawn from the 
three different Monte-Carlo samples are calculated. In order 
to do so, the distribution of mmax around each observed 
cluster mass, /x(mmax : Med), is examined for the three 
samplings, whereby we only use those that have a max- 
imal stellar mass, mmax* = 150 M©. If the mean value, 
?Timax(M(,ci), of the Montc-Carlo distribution, /i(mmax), for 
a cluster mass is higher than the observed maximal mass, 
JTimaxobs, thcu the distribution is integrated from the lower 
limit, miow, to mmaxobs and divided by the integral from 

f n ^TTT^^^" sort 
-P(mniax ^ mniaxobs) = .—Tin^con sort • (24) 

/ u dm 

This is the probability of observing a maximum mass mmax 
as small as or smaller than mmaxobs- 

If, on the other hand, the mean value is smaller than 
mmaxobs, the integral is taken from mmaxobs to the upper 
mass limit, mmax*, and is divided by the integral from the 
mean value, 'mmai''""' to mmax*, 

-P(mmax ^ mmaxobs) ~ fTTimax* j ' (2^) 

This is the probability of observing an mmax as large as 
or larger than mmaxobs- Together, eas. 1241 and 1251 are the 
probability of observing an mmaxobs such that | mmaxobs — 
THmax (Mecl)l ^ I (Mcci)|. Figs.irniandlTTlshow 

a schematic view of this integration process by means of two 
examples (Taurus- Auriga and NGC2264) for all three cases 
of sampling. 

The resulting combined probabilities (11 in Tab. |5| ac- 
quired by multiplying the individual probabilities) for all 
clusters are 10~® for random sampling, 10~® for constrained 
sampling and 10~* for sorted sampling, with the highest 
probability thus being attained for sorted sampling. Tab.|5| 
shows the individual probabilities P for each mmaxobs- 

To obtain a completely different statistic on the cor- 
respondence between dat a and theory, we also apply a 
Wilcoxon signed-rank test jBhattacharvva fc Johnsonll977ll 
to determine how significant the differences between the 
data and the 'mmai"""' """^^ relations are. For this test the 
differences of the data points and the mean values are cal- 
culated, 

A I ran. con, sort / 1- ■ \ /or-x 

mj = mmaxobs, j - mj^ax,j (tOr a glVCU Med, j ) , (26) 

and then ranked according to their absolute value. Af- 
terwards, only the positive-signed ranks are added and 
this sum is then cross-checke d with tabulated values 
ijBhattacharvva fc ,Tohnsonlll977t) in order to get the prob- 
ability P that the null hypothesis (data points and the 



Table 2. Probabilities that the rrimaxobs from one of the 
three Monte-Carlo samples. 



Cluster 


random 


constrained 


sorted 




sampling 


sampling 


sampling 




% 


% 


% 


Tau-Aur 


25.0 


36.2 


45.4 


Ser SVS2 


12.5 


33.1 


41.5 


IN VjV_j J.OOO 


24.7 


29.4 


41.1 


p Oph 


31.5 


60.6 


81.7 


IC348 


9.2 


38.2 


52.5 


NGC2024 


69.1 


86.0 


76.4 


o- Ori 


69.1 


86.0 


76.4 


Mon R2 


25.8 


27.1 


37.5 


NGC2264 


57.6 


81.9 


84.7 


NGC6530 


27.2 


32.4 


44.3 


Ber 86 


39.3 


53.9 


76.3 


M42 


43.9 


53.6 


75.5 


NGC2244 


34.1 


40.4 


68.0 


NGC6611 


6.4 


5.6 


32.2 


Tr 14/16 


42.4 


44.8 


96.6 


Arches 


75.0 


88.2 


67.4 


R136 


65.7 


82.3 


42.2 


n 


3.37 X 10-9 


9.22 X lO"'^ 


9.18 X lO"'^ 



''Timax'^""' relations are the same within the uncertain- 
ties) is true. In the case of random sampling (§ 12. 2. in and 
constrained sampling (§ ITT^ P = 0.00014, while P = 
0.0458 for sorted sampling (^ HTHi . 

Thus, both tests taken together suggest strongly that 
sorted-sampling best describes the empirical data. The phys- 
ical interpretation of this result is discussed in § 2] 



3 THE MONTE-CARLO SIMULATIONS OF 
THE IGIMF 

We have thus seen that the observational data strongly 
favour a particular mmax(Afed) relation, namely the 
mmax ~ '™-max relation. This has profound implications 
for composite stellar populations. 

Fig. I12 | shows the r esult of the semi-analytic approach 
by iKroupa fc Weidneil l)2003h assuming (3 = 2.35 for the 
ECMF. The stellar IMF in each cluster has, in all cases, the 
standard or canonical three-part power-law form (eq. |5J. 
For the minimum "cluster" mass, Mcci,min = 5 M© (a 
dozen stars), and for the maximal cluster mass, Med, max = 
IO'^Mq, are used. The power-law index qigimf of the 
resulting semi-analytical IGIMF is well approxima t ed by 
QiGiMF = 3.00 for m> 1M0. In lKroupa fc Weidneil tOOm 
we already noted that this is p robab l y the reason why the 
Galactic-field-IMF deduced bv IScald Jl986D (as « 2.7) is 
steeper than the canonical stellar IMF. 

We now apply the Monte-Carlo experiments introduced 
above to the computation of the IGIMF (eq. The result- 
ing IGIMF is constructed by mass-binning all stars in all 
clusters. It is shown as a long dashed line in Fig. ll3l for con- 
strained sampling and as a short dashed line for sorted sam- 
pling. Additionally the 'input' standard stellar IMF (solid 
line) and the semi-analytical IGIMF from Fig. 1121 [dotted 
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log 

10 

Figure 10. All panels: Vertical slice through Figs. ISll^andl^at 

Meei = 25 Mq (Taurus- Auriga) . The vertical axis plots the num- 
ber of clusters /i which have the maximal mass rrtmax in a logarith- 
mic bin AlogjQmmax = 0.01044. Top panel: Histogram of the 
number of clusters for the random samphng models (S I2.2.1I . The 
vertical dashed line marks the observed upper mass in Taurus- 
Auriga, while the dotted line is the mean (or expectation) value 
of the Monte-Carlo experiment, 'm:™^(M(,ci). The probabihty 
of observing small as or smaller than ^rij^axobs is 

^'("^max ^ '"maxobs) which is the area under the curve left of 
the dashed line divided by the area under the curve left of the 
dotted line. Middle panel: Like top panel but for the constrained- 
sampling Monte-Carlo experiment (§ 12.2.21 . Bottom panel: As 
top panel but in this case for sorted samphng (^ l'2.2.',H .All panels: 
Note that the areas to the left and right of the mean value do not 
appear equal because the binning was done linearly but is plotted 
on a logarithmic scale. 



line) axe shown. The IGIMF obtained with sorted samphng 
agrees very well with the semi-analytical result, while the 
IGIMF obtained from constrained sampling shows a flatter 
slope. Even this IGIMF nevertheless differs strongly from 
the standard IMF with a Salpeter slope, as a result of the 
imposed condition of a given cluster mass. 



3.1 Different ECMF 



Follow ing- on from our discussion in IWeidner fc Kroupal 
J2005I) we explore how a different ECMF affects the IGIMF. 
Any ECMF with /3 > 2.35 will increase the steepening of 
the IGIMF. However, below about 50 or 100 Mr the E CMF 
is poorly defined observationally jLada fc Ladall2'003r) . and 
we consider here the implication of a flatter ECMF at low 
masses while keeping /3 = 2.35 for Med > 50 or 100 Mq. 




logio(mn,ax/Mo) 

Figure 11. All three panels like Fig. 1101 but in this case for a 
cluster mass of 355 Mq (NGC2264). For the sorted sampling 
{bottom panel) the mean is lower than the observed value of mmax- 
Therefore the probability is calculated by dividing the area to the 
right of the observed value {dashed line) through the area to the 
right of the expectation value {dotted line). 



In Fig.^]we explore the effect of a different ECMF than 
our standard assumption (/3 = 2.35 over the whole range of 
cluster masses) . This flgure is the same as Fig. 1131 but the 
constrained sampling and sorted sampling results for two 
EMCFs with /3i = 1 for clusters below 50 M© and 100 Mq, 
and (32 — 2.35 for clusters above these masses, are shown as 
dotted and dashed lines, respectively. The steeper ones show 
the sorted-sampling cases. The change of the ECMF reduces 
the effect of the ECMF on the IGIMF but is still clearly 
visible, giving aiGiMF,con = 2.56, aiciMF.sort = 2.88 for /? = 1 
(Mod < 50 Mq) and qigimf.coii = 2.52, aiGiMP.sort = 2.86 
for /3 = 1 for (Med < 100 A/q), as opposed to our standard 
assumption for which aioiMF.con = 2.8 and aiciMF.sort = 
3.0. 



4 DISCUSSION AND CONCLUSIONS 

In this contribution a number of Monte-Carlo experiments 
are conducted in order to constrain the relation between the 
maximal mass a star can have in a cluster and the mass 
of the cluster, and to further study the IMF of composite 
populations, ie. the integrated galaxial IMF. 

We consider three possible ways of forming clusters: 
(1) Completely randomly. Clusters are fllled randomly with 
stars and then their masses, Med, are calculated (random 
sampling). (2) Cluster-masses are picked from an ECMF 
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Figure 12. Solid line: Canonical stellar IMF, ^(m,), in logarith- 
mic units and given by the standard three-part power-law form 
(eq. |5J, which has 03 = 2.35 for m > 0.5 Mq. Dotted line: semi- 
analytical ^igimfC™) for /3 = 2.35. The IMFs are scaled to have 
the same number of objects in the mass interval 0.01 — 1.0 Mq. 
Note the turn down near mmax* = 150 Mq which comes from 
taking the fundamental up per mass limit explicitly into account 
iWeidner fc Kroupal [2004) . Two lines with slopes anac = 2.35 
and aiine = 3.00 are indicated. 
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Figure 13. Solid line: Standard stellar IMF with ag = 2.35 for 
m > 1 M© (same as in Fig. Dotted line: IGIMF resulting 

from the semi-analytical approach with /3 = 2.35 (as in Fig. 1121 . 
Short dashed line: IGIMF obtained from sorted sampling. Long 
dashed line: IGIMF produced by constrained sampling of stars. 
As the IMF below 1 Mq does not change, only the region above 
1 Mq is plotted here. 
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Figure 14. As Fig. 1131 fall lines from Fig. 1131 are solid) but in 
addition the results for two different ECMFs are plotted as dotted 
and dashed lines. In the dotted case the slope of the ECMF is f3i = 
1 for clusters below 50 Mq . In the dashed case /3i = 1 below 
100 Afg. In both cases P2 = 2.35 for clusters above these limits. 
As the IMF below 1 Mq does not change, only the region above 
1 Mq is shown. The respective models assuming sorted sampling 
are always steeper than those assuming constrained sampling. 



and used as a constrain in constructing the stellar content 
of each cluster (constrained sampling). (3) Cluster-masses 
are picked from an ECMF, and the clusters are then filled 
with stars by randomly selecting from the canonical IMF, 
sorting the stellar masses in ascending order and constrain- 
ing their sum to be the cluster mass (sorted sampling). In 
all cases (1-3), the most massive star, rrimax, in each clus- 
ter is found, and the average or expectation value, mmax, is 
calculated for the ensemble of clusters near Mod to give the 
relations "m™i'=°"' (Meci) , where "1, u" refers to models 
with or without a physical stellar mass limit of 150 Mq . 

The most important and surprising result is that the 
sorted-sampling algorithm best represents the observational 
data of young (^ 3 Myr) clusters. Constrained and random 
sampling do not fit the observations. 

That our sorted-sampling algorithm for making stars 
fits the observational maximal-stellar-mass — star-cluster- 
mass data so well would appear to imply that clusters form 
in an organised fashion. The physical interpretation of the 
algorithm (i.e. of the Monte-Carlo integration) is that as a 
pre-cluster core contracts under self gravity the gas den- 
sities increases and local density fluctuations in the tur- 
bulent medium lead to low-mass star formation, perhaps 
similar to what is seen in Taurus-Auriga. As the contrac- 
tion proceeds and before feedback from young stars begins 
to disrupt the cloud, star-formation activity increases in 
further density fluctuations with larger amplitudes thereby 
forming more massive stars. The process stops when the 
most massive stars that have just forme d supply sufficien t 
feedback energy to disrupt the cloud (lElmegreenI Il983f) . 
Thus, less-massive pre-cluster cloud-cores would "die" at 



12 C. Weidner and P. Kroupa 



a lower maximum stellar mass than more massive cores. 
But in all cases stellar masses are limited by the physical 
maximum mass, m ^ mmax(Meci) ^ mmax*. This sce- 
nario is nicely consistent with the hydrodynamic cluster 
forma tion calculations presented by [ Bonnell. Bate fc Vina 
J2003D andlBonnell. Vine fc Batd J2004) . We note here that 
lBonnSL~V^i^^Bat^ ~ il20o3^ their theoretical clus- 

ters to form hierarchically from smaller sub-clusters, and 
together with continued competitive accretion this leads to 
the relation mmax oc M^^^ (eq. 1131 in excellent agreement 
with our compilation of observational data for clusters with 
masses below M^d ~ 4000 Mq . While this agreement is 
stunning, the detailed outcome of the currently available 
SPIf modelling in terms of stellar multiplicities is not right 
jGoodwin fc Kroupal lioosll . and feedback that ultimately 
dominates the process of star- formation, given the generally 
low star- formation efficienci es observed in cluster-forming 
volumes jLada fc Ladall2003D . is not yet incorporated in the 
modelling. 

Stellar evolution is the major caveat here. But the com- 
parison of different models (see Appendix|nj shows a general 
agreement of the lifetimes and relevant parameters (mass, 
Tcff and luminosity) for the models considered here. There- 
fore, not different models but an intrinsically steeper IMF 
{as 2.8) could shift the expectation values for random and 
constrained sampling into the observed regime. Such a steep 
IMF may be possible if it is masked by unresolved multiple 
stars, something we are investigating now. 

Furthermore, the Monte-Carlo expe ri ment s as- 
certain the result s of i Kroupa fc WeidneJ ll2003f) and 
IWeidner fc Kroupal J2005I) regarding the steep IGIMF, 
especially so if sorted sampling is used. In the constrained 
sampling case the IGIMF slope is still steeper than the 
input IMF but less steep than with sorted sampling. 
But it should be noted h ere that a very recent result by 
lElmegreen fc Scald (l2006l) shows that it is also possible to 
interpret PDMF variations falsely, as IMF variations when 
the SFR is assumed to be constant when in reality being 
burst-like. This result has not yet been implemented in our 
description of the IGIMF. 

In summary: 



- There exists a well-defined relation, mmax = 
mmax(Afeci), between the most-massive star in a cluster and 
the cluster mass. The conjecture that a cluster consists of 
stars randomly picked from an invariant IMF between 0.01 
and 150 Mq would therefore appear to be wrong. 

- Star clusters appear to form in an ordered fashion, 
starting with the lowest-mass stars until feedback is able 
to outweigh the gravitationally-induced formation process. 

- IGIMFs must always be steeper for m > 1 Mq than the 
stellar IMF that results from a local star-formation event. 



It will be important to further test the results presented 
here on the mmax(A/oci) relation by compiling larger ob- 
servational samples of young clusters. As this contribution 
has shown, it appears that the mmax(Afeci) relation would 
be rather fundamental to galactic and extragalactic astro- 
physics. 
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Table Al. Observed mass ranges and the corresponding mean 
masses for those clusters for which we need to derive cluster 
masses. 



Cluster 


number 


% of 


total 


mass range 








of stars 


all 


number 


[Me] 


[Me] 


[Mq] 




observed 


stars 


of stars 








Tau-Aur 


100 


92.8 


108 


0.02 - 2.2 


0.236 


25.5 


IC348 


241 


62.5 


386 


0.08 - 6.0 


0.282 


109 


Mon R2 


475 


55.2 


861 


0.1 - 10.0 


0.301 


259 


NGC2264 


600 


55.3 


1086 


0.1 - 25.0 


0.327 


355 


Ber 86 


340 


7.7 


4421 


0.8 - 40.0 


0.338 


1500 


NGC2244 


400 


2.2 


17937 


2.0 - 70.0 


0.348 


6240 


NGC6611 


362 


0.6 


56563 


5.0 - 85.0 


0.352 


20000 


Tr 14/16 


768 


0.6 


120000 


5.0 - 120.0 


0.357 


43000 
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New Astron., 7, 395 
Zhang Q., Fall S. M., 1999, ApJL, 527, L81 
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by the author. 



APPENDIX A: CLUSTER AND MAXIMAL 
STAR MASS DETERMINATION 

The masses of the clusters in Tab. Q are acquired as follows: 
In the case of NGC1333, NGC2024, NGC6530, M42, p 
Ophiuchi, a Orionis, Serpens SVS2, Arches and R136 the au- 
thors of the corresponding papers provide the required mass 
estimates. The masses for NGC2244, NGC2264, NGC6611, 
IC348, Monoceros R2, Taurus- Auriga, Berkley 86 and Trum- 
pler 14/16 are calculated by determining the fraction (given 
as a percentage in Tab. lAH of observed stars in compari- 
son with a canonical IMF from 0.01 Mq up to the observed 
upper mass limit. With the fraction of all stars the total 
number of stars in the cluster is estimated by dividing the 
observed number of stars by the fraction. The total mass of 
the cluster, Meci, is then calculated by multiplying the total 
number of stars with the mean stellar mass, mmean, for the 
canonical IMF from 0.01 Mq to the observed upper mass 
limit. The relevant values are shown in Tab. lAll 

For the maximal stellar masses in these clusters the val- 
ues within the papers are used whenever possible, which 
is for NGC2244, NGC2264, NGC6611, M42, p Ophiuchi, 
a Orionis, Monoceros R2, Berkley 86, Trumpler 14/16, 
Arches and R136. In the other cases (NGC1333, NGC2024, 
NGC6530, IC348, Serpens SVS2, Taurus-Auriga) mass es- 
timates are derived from the spectral types of the most lu- 
minous me mbers using the spectral-type mass-relation from 
ICoxU2000h . 



APPENDIX B: FITTING FORMULAE FOR 
MASSIVE STAR EVOLUTION 

Because the lHurlev. Pols & Touli (l2000ll single stellar evolu- 
tion (SSE) package is only calibrated for stellar models up 
to 50 Mq , additional fitting formulae have been developed 
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for more-massive stars. Based on the lSchaller et aTl (Il992f) 

models for 60, 85 and 120 Mq, functions for m{t), L(t) and 
Tcs(t) have been obtained: 



Bl Mass Evolution 

As long as the age, t (in Myr), of the star is below Tm, 
the main-sequence life-time, the mass-evolution can be de- 
scribed according to, 



m(t) — mini ■ e '"i'*' . 

For Tm < t ^ Tn^ + dt 

at 

For both the parameters are 

' ■tam\ 

+ tab 



ai 

h 

bl 



(Bl) 
(B2) 



106 



= mini • (1 - /i), 



{frr 

02 ■ ( 1 + 



ft, 



^ dt a2 ) 



and constants 



rriaa 


= 48.0, 




rriab 


= 24.7, 




rriac 


= 3.15, 




dt 


= 0.42, 




fm 


= 3.523808 • 


10 


ft 


= 6.190428 ■ 


10 


taf 


= 23.5, 




tam 


= 1.25 • 10^ 




tab 


= 2.5 • 10^ 





When t is larger than Tm + dt the star is considered dead. 
No remnant mass, Toff or luminosity is assigned. 

The resulting curves for a 120, 85, 60 and 50 Mq star 
{solid lines) in comparison with the model data {dotted lines) 
are plotted in Fig. IBll 



B2 Effective Temperature Evolution 

For t ^ Tm the equation 



120 




2 3 
age/Myr 

Figure Bl. Stellar mass ev olution for massive s tars. The dotted 
lines are the Geneva models iSchaller et al.ll992^ . while the solid 
lines arc the results from the fitting formulae described here. A 
ISurlcv. Pols Tout. L2noa'> 50 M0-star is shown as a dashed line. 



as = 



/2 = 



ttai mini 5S 85M0, 

tta2a ■ mini + tta2h mini < 85M0, 

(^peak '?low) 



dt 



Tlow — ^?cfT, in 

and the constants 



I -ISL 



Tim-^ 

Tib, 

71ni2 
Tih2 

tta\ 
tta2h 

^pcak 



4.3 • 10"^ 
4.425, 

7.3333 • 10~*, 

4.64, 

2.3, 

9.0, 

-0.182857142, 

24.54285719, 

4.8. 



The resulting Tcff fitting curves are plotted in Fig. 
for the same masses as in Fig. IBll 



ns{t) ^Teff.ini-el , 

adequately captures the evolution. 

For Tm < t ^ n-n + dt, 



ns{t)= f2-t+{Ti,„-{f2 

is used. 



0), 



The parameters are 

Tcff.ini ~ 



(B3) 



(B4) 



Timi • mini + Tibi mini ^ 6OM0, 
Tima • mini + Tiba mini < 6OM0 , 



B3 Luminosity Evolution 

The luminosity evolution is divided into three parts. For 

t 5^ Tm , 



L{t) = ( 



-^jump Lini 



For Tm < t < Th 



■ t + Lini. 



(B5) 



L{t) = I/jump + (Lpcak - -^jump) ' sin(Lsin ' {t - Tm)) ■ (B6) 

And finally for Th ^ t ^ {Tm + dt) 
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age/Myr 

Figure B2. Effective temperature evolution for massive stars. 
Tlie line styles are as in Fig. lBll As the Tcff for massive stars are 
rather similar, the lines in the plot have been shifted upwards as 
indicated in the plot. 
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Figure B3. Luminosity evolution for massive stars. The line 
styles are as in Fig. IBl| 



Fig. IB3I shows the results of the luminosity- fitting for- 
mulae. 



L{t) = ai-t + b^. 

Here the parameters are 



(B7) 



Ml 

-^'jump 
-^'pcak 

63 



mini 

Ml ' 



exp 



1 • 10~ 



Tbb 



min, > 



-^jump 



^low 



Tm -\- dt ~ Tb 
L\ow - (^^4 ■ (Tm + dt)). 



The constants are 




= 0.145590532, 


mlh 


= 1.722994092, 


Ljc 


= 9.0 ■ 10"^ 


Lpi 


= 2.5, 


Ljc 


= 3.95, 


Tbb 


= 4.6296293 • 10 


Tbm 


= 1.1111 ■ 10"^ 




= 15.7, 




= 1.35, 


Lpc 


= 4.215, 




= 5.2. 



APPENDIX C: COMPARISON OF DIFFERENT 
MODELS 

Several different sets of theoretical models for stellar evo- 
lution of massive stars exist. Figs. lUlT IC2I and IC3I compare 
the mass, T^s and luminosity evolution of three different 
sets of models ("Schaller et al.ll992l : lMevnet fc Maederl2003l : 
[Hurlcv. Pols & Tout 200^ . These models agree qualitatively 
on the compared stellar properties but show minor differ- 
ences in the details. 
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age/Myr 

Figure CI. Stellar mass evolution fo r massive stars. Th e 
solid lines are the Geneva models iSchaller et al.l Il992t) . 
while the das hed lines arc the Geneva models with rota- 
tion iMevnet fc Maeder 2003, thick dashed line: no rotation, 
medium dashed line: 300 km/s, thin dashed line: 500 km/s, 
only for m = 60 Mp ) an d the results from the SSE package 
jHurlev. Pols fc Toujl200oh are shown as dotted lines. 



O 



5.5 




2 3 
age/Myr 

Figure C3. Luminosity evolution for massive stars. The line 
styles arc as in Fig. IGll 




age/Myr 



Figure C2. Effective temperature evolution for massive stars. 
The line styles are as in Fig. IGll As the T^g for massive stars are 
rather similar, the lines in the plot have been shifted (as indicated 
in the plot) in the following way: the upper-most by -1-1.0 dex and 
the second one by -1-0.5 dex. The lowest plot has not been shifted. 



